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Abstract. We report about a recent calculation of the heavy quarkonium mass and decay width in 
a quark-gluon plasma, whose temperature T and screening mass mo satisfy the hierarchy ma s ^> 
%T ^> ma s 2 ^> mo, m being the heavy-quark mass, up to order ma s 5 . The calculation may be relevant 
to understand the behavior of the T(15) in a quark-gluon plasma at present-day colliders. 
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INTRODUCTION 

The behavior of quarkonium in a medium has been the object of many studies since 
it was realized that the medium may induce quarkonium dissociation and that the 
dissociation pattern of the different quarkonia may provide an effective thermometer 
of the medium [1]. How the medium influences the quarkonium and eventually leads to 
its dissociation is the problem that such studies address. 

In recent years, new insight in the problem has been gained by studying low-lying 
quarkonia in a weakly-coupled quark-gluon plasma. Low-lying quarkonia, like the bot- 
tomonium and charmonium ground states, are (believed to be) characterized by a typical 
inverse radius, which is of the order of moc s , larger than the non-perturbative scale of 
QCD, Aqcd> and a typical energy E of the quarks in the bound state, which is of the 
order of ma s 2 , larger than or of the same order as Aqcd- Such states are Coulombic. In 
the case of the bottomonium ground state, it is often also assumed that E 3> Aqcd- In a 
weakly-coupled plasma, both the temperature T of the plasma and the screening mass, 
mo, are larger than Aqcd- Moreover, the perturbative power counting implies that mo 
is suppressed with respect to T, for m D shows up at one loop. With these specifications, 
observables describing low-lying quarkonia in a weakly-coupled quark-gluon plasma 
may be calculated to a large extent analytically in perturbation theory, which makes the 
interest of these systems. As we will argue later, this could be the case of the T(15) at 
the LHC. 

In [2, 3, 4, 5], the perturbative quarkonium static potential was studied for quark- 
antiquark distances r such that T ^> 1/r > mo- The surprising output of these studies 
has been that the imaginary part of the gluon self energy, which may be traced back 
to the Landau-damping phenomenon, induces an imaginary part of the static potential 
and hence a thermal width of the static quark-antiquark bound state. In [6], static 
particles in a hot QED plasma were considered in the situation 1/r ~ mo, confirming 



previous results. In [7], the quarkonium static potential was first studied in an effective 
field theory (EFT) framework that exploits systematically the hierarchy of different 
energy scales in the problem. The static potential was studied for distances that range 
from larger to smaller than \/T and a new dissociation mechanism was identified. 
This is the color-singlet to color-octet break-up mechanism that provides the dominant 
contribution to the thermal decay width when E 3> mo- In [8], moreover, the relation 
between the quarkonium potential and the correlation of two Polyakov loops, a quantity 
often evaluated on the lattice, has been investigated. A comprehensive study of non- 
relativistic bound states in a hot QED plasma in a non-relativistic EFT framework was 
performed in [9, 10]. It has been pointed out that quarkonium melts at temperatures 
larger than the melting temperature, r me i t i ng , defined as the temperature for which the 
thermal decay width becomes as large as the binding energy; parametrically, it holds 
that 7rr me i t j ng ~ mg 4//3 . Since this temperature is lower than the temperature at which 
screening sets in, i.e. the temperature at which mo is of the size of the inverse of 
the radius of the bound state, this finding challenges the long-time accepted view that 
quarkonium remains dissociated in a medium as long as color is screened; in fact, the 
existence of the melting temperature implies that quarkonium remains dissociated at 
lower temperatures, experiencing a larger dissociation rate. 

In the following, we will report about the results of [11], where, in a specific range 
of temperatures, the spectrum and widths of quarkonia up to order ma s 5 have been 
computed. In the assumed range of parameters, the scales of the bound state (m, ma s 
and ma s 2 ) and the thermodynamical scales (T and mo) fulfill a specific hierarchy that we 
will exploit either by constructing a corresponding hierarchy of EFTs or by subsequently 
integrating out different momentum regions. 

SCALES AND EFTS 

We consider quarkonium in a thermal bath, whose energy scales satisfy the following 
hierarchy: 

m > ma s > KT > ma s 2 > m D . (1) 

This implies that mg 3 >T> mg 4 and that nT is lower than 7rr me i t i ng , i.e. quarkonium 
exists in the plasma. We will further assume that all these scales are larger than Aqcd and 
that a weak-coupling expansion is possible for all of them. Finally, in order to produce 
an expression for the spectrum that is accurate up to order ma s 5 , we will assume that 
[m D /(ma s 2 )] 4 <g. 

It is not clear if this hierarchy is realized by any of the quarkonium states at present 
heavy-ion facilities. Possible candidates are the bottomonium 15 states (T(15), rj/,) at 
the LHC, for which it may hold m b ps 5 GeV > m h a s « 1.5 GeV > nT w 1 GeV > 
ma s 2 « 0.5 GeV > m D . 

It is useful to work in the so-called real-time formalism, which amounts to modifying 
the temporal integration in the partition function to include real times. As a consequence, 
the degrees of freedom double. However, the doubling of the degrees of freedom affects 
only gluon loops (at the order we are working only the gluon self energy), while physical 
heavy quarks decouple from unphysical ones [7], so the price to pay is minimal. The 



advantage is that the real-time formalism allows a treatment of the quarkonium in the 
thermal bath very similar to the EFT framework developed for zero temperature [12]. 
Many results obtained there may, indeed, be simply translated here. 

The EFT that follows from QCD by integrating out gluons and quarks of energy or 
momentum of order m in the quark-antiquark sector is non-relativistic QCD (NRQCD) 
[13]. The EFT that follows from NRQCD by integrating out gluons of energy or mo- 
mentum of order ma s is potential non-relativistic QCD (pNRQCD) [14]. According to 
the hierarchy (1), both these scales are larger than T, which, therefore, may be set to 
zero in the matching to the EFT. As a consequence, the Lagrangians of NRQCD and 
pNRQCD are the same as at zero temperature. 

Integrating out T from pNRQCD modifies pNRQCD into hard-thermal loop (HTL) 
pNRQCD, pNRQCD HT L, U, 15]. With respect to pNRQCD, the pNRQCD HT L La- 
grangian gets relevant modifications in two parts. First, the Yang-Mills Lagrangian gets 
an additional HTL part [16]. This, for instance, modifies the longitudinal gluon propa- 
gator in Coulomb gauge into (k 2 = k 2 ) 
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where "+" identifies the retarded and "— " the advanced propagator. Second, the po- 
tentials get in addition to the Coulomb potential, which is the potential inherited from 
pNRQCD, a thermal part, 8V. 



POTENTIAL, ENERGY AND DECAY WIDTH 

In the following, we will provide the thermal corrections to the color-singlet quark- 
antiquark potential, the thermal corrections to the spectrum and the thermal decay width, 
aiming at a precision of the order of ma s 5 . 




FIGURE 1. The single line stands for a quark-antiquark color-singlet propagator, the double line for a 
quark-antiquark color-octet propagator and the circle with a cross for a chromoelectric dipole interaction. 



Integrating out the scale T 

The relevant diagram contributing to the potential is shown in Fig. 1. It reads 

where D^ v stands for the gluon propagator, h () = p 2 /m+ 0C S / (6r) is the octet Hamilto- 
nian and the loop integral has been regularized in dimensional regularization (D = 4 + £ 



and fi is the subtraction point). In the loop integral, we integrate over the momentum 
region k Q ~T and k^T. Since T 3> (E — h ), we may expand 
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FIGURE 2. Detail of the different terms contributing to the potential. 
The real part of the thermal correction to the color- singlet potential reads 



ReSV v (r) = 



— a 2 rT 2 A a s T 2 

9 s 9m 



4a s / r 



9a s 3 17 a 2 



+ ■ 



4 7ra 



37T 8 r 3 mr 2 9 m 

-2C(3)^r 2 rm 2 +|c(3)« 2 r 2 r 3 , 
7T 3 

IT 2 5 
/ r = -+ln-2 -7 £ + ln(4^)--, 
e ju z 3 

W 2 ) 



fi'(r) + 



? 



(4) 
(5) 
(6) 



where the first line of Eq. (4), which is of order g 2 r 2 T 3 xE/T, comes from diagram a) 
in Fig. 2, the second one, which is of order g 2 r 2 T 3 x (E/T) 3 , comes from diagram b), 
and the third one, which is of order g 2 r 2 T 3 x (mo/T) 2 , comes from diagram c); n/ is 
the number of light quarks. 

The imaginary part of the color- singlet potential, which comes from diagram c) in 
Fig. 2, reads 

2 o o / 2 T 2 2 C'(2)\ 
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This contribution, which may be traced back to the Landau-damping phenomenon, is of 
order g 2 r 2 T 3 x (mo/T) 2 . 

Evaluating Re 8V s (r) and Im 8V s (r) on a quarkonium state with quantum numbers 

(T) 

n and /, we obtain the thermal correction to the energy, 8E K n / , and the thermal width, 
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Integrating out the scale E 

The diagram shown in Fig. 1 also carries contributions coming from the energy scale 
E. They may be best evaluated in pNRQCDhtl by integrating over the momentum 
region ko ~ E and k ~ E and using HTL gluon propagators. Since k ~ E <C T, we 
may expand the Bose-Einstein distribution 

nB(k) = l-- l - + -L + ...; (10) 

moreover, since k ~ £ ^> mo, the HTL propagators can be expanded in m 2 D jE 2 <C 1. 

The momentum region ko ^ E and A: ~ £ is characterized by two possible momentum 
sub-regions. This can be understood by considering the integral 
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For ko ^ E and k ~ E, it exhibits an off-shell sub-region, ko — k^E, and a collinear 
sub-region, ko — k ~ m 2 D /E. Note that, according to (1), the collinear scale satisfies 
mg A 3> mjy/E 3> mg 6 , i.e. it is smaller than mo by a factor of mo/E <C 1 but still larger 
than the non-perturbative scale g 2 T by a factor T/E ^> 1. 

The thermal correction to the energy, 8E^, coming from the scale £, reads 

5£ lf = -^a s rm 2 a^ 2 [5n 2 + l-3/(/ + l)]. (12) 

We note the complete cancellation of the vacuum contribution (which includes the Bethe 
logarithm) against the thermal contribution originating from the "—1/2" term in the 
expansion of the Bose-Einstein distribution (see Eq. (10)). 




The thermal width, T y n coming from the scale E, reads 
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where /^o = —0.49673, /2,o = 0.64070, .... The leading contribution is given by the first 
three terms, which are of order CC^T. This contribution to the thermal width is generated 
by the possible break up of a quark-antiquark color-singlet state into an unbound quark- 
antiquark color-octet state: a process that is kinematically allowed only in a medium. The 
singlet to octet break up is a different phenomenon with respect to the Landau damping. 
In the situation E ^> m D , the first dominates over the second by a factor (ma 2 /m D ) 2 . 



The diagram shown in Fig. 1 also carries contributions coming from the energy scale 
mj). These contributions are suppressed with respect to the other terms calculated. 



The thermal corrections to the spectrum and the thermal decay width develop diver- 
gences at the different energy scales. These are artifacts of the scale separations and 
cancel in the final (physical) results. 

Concerning the thermal decay width, the divergence at the scale ma 2 in Eq. (13), 
which is of ultraviolet (UV) origin, cancels against the infrared (IR) divergence at the 
scale T in Eq. (9). 

In the spectrum, the pattern of divergences is more complicated and is summarized 
in Tab. 1 . The table may be read vertically or horizontally. If read vertically, it shows 
a typical EFT cancellation mechanism: at the scale moc s we have non-thermal IR diver- 
gences in the potentials, these cancel against non-thermal UV divergences at the scale 
ma 2 [17], the non-thermal contribution at the scale T is scaleless and vanishes in dimen- 
sional regularization; thermal IR divergences at the scale T cancel against thermal UV 
divergences at the scale ma 2 . If read horizontally, it shows a cancellation mechanism 
that is familiar in thermal field theory: at the scale T, thermal IR divergences cancel 
against non-thermal IR divergences while non-thermal UV divergences cancel against 
IR divergences that appear in the potentials at the scale ma s ; at the scale ma 2 UV ther- 
mal divergences cancel against UV non-thermal divergences. The cancellation between 



Integrating out the scale mo 



CANCELLATION OF DIVERGENCES 



TABLE 1. The pattern of IR and UV diver- 
gences in the quarkonium spectrum at different 
energy scales. The final result is finite. 
Scale I Vacuum i Thermal 
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non-thermal and thermal parts is possible because the latter may carry temperature in- 
dependent terms (see, for instance, the "—1/2" term in Eq. (10)). Note that both at the 
scales T and ma 2 , the spectrum is finite. 



SUMMARY 



For a quarkonium state that satisfies the hierarchy specified in Eq. (1) and in the follow- 
ing discussion, the complete thermal contribution to the spectrum up to &{mal) is ob- 
tained by summing Eqs. (8) and (12) and subtracting from the latter the zero-temperature 
part, this gives 



8E 



(thermal) 



'n.l 



— a 2 T 2 a 



3n 2 -/(/ + !) + 



-2y E 



32 5,, 



1 



' 27 n 3 n(2/ + l) 



16 1 27 



128£„a s 3 



SlTt 

+aln 2 [5n 2 +l-3/(/+l)j 



>> + t 



a s Tm 2 D 



+ ^C(3)« s 2 r 3 }, 



(14) 



where L n j are the QCD Bethe logarithms: L lfi = -81.5379, L 2 ,o = -37.6710, . . . [18]. 

For a quarkonium state that satisfies the hierarchy specified in Eq. (1) and in the 
following discussion, the complete thermal width up to ^(ma s 5 ) is obtained by summing 



Eqs. (9) and (13), this gives 
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As a qualitative summary, we observe that, at leading order, the quarkonium masses 
increase quadratically with T, which implies the same functional increase in the energy 
of the leptons and photons produced in the electromagnetic decays. Electromagnetic de- 
cays occur at short distances ~ 1/m <C l/T, hence the standard NRQCD factorization 
formulas hold. At leading order, all the temperature dependence is encoded in the wave 
function at the origin. The leading temperature correction to it can be read from the po- 
tential and is of order n 4 T 2 / (m 2 a s ). Hence, a quadratic dependence on the temperature 
should be observed in the frequency of produced leptons or photons. Finally, at leading 
order, a decay width linear with temperature is developed. The mechanism underlying 
this decay width is the color-singlet to color-octet thermal break-up, which implies a 
tendency of the quarkonium to decay into a continuum of color-octet states. 
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